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ABSTRACT 


The emissivity for the neutrino pair synchrotron radiation in strong magnetic 
fields has been calculated both analytically and numerically for high densities 
and moderate temperatures, as can be found in neutron stars. Under these 
conditions, the electrons are relativistic and degenerate. We give here our 
results in terms of an universal function of a single variable. For two different 
regimes of the electron gas we present a simplified calculation and compare our 
results to those of Kaminker et al. Agreement is found for the classical region, 
where many Landau levels contribute to the emissivity , but some differences 
arise in the quantum regime. One finds that the emissivity for neutrino pair 
synchrotron radiation is competitive, and can be dominant, with other neutrino 
processes for magnetic fields of the order B rsj 10 14 — 10 15 G. This indicates the 
relevance of this process for some astrophysical scenarios, such as neutron stars 
and supernovae. 


Subject headings : Stars : Magnetic Fields — Stars: Neutron. 
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1. INTRODUCTION 


Estimates of the magnetic field strength at the surface of neutron stars are obtained 
from several different scenarios: theoretical models of pulsar emission ([Ruderman 1972 ), the 
accretion flow in binary X-ray sources ([Ghosh and Lamb 19781 ) and observation of features 
in the spectra of pulsating X-ray sources which have been interpreted as cyclotron lines 
( [Triimpcr ct al. 1978| , IWheaton et ah 1979| , IGruber et al. 19801 , |Mihara ct al. 1990|) . For 
a recent review of all these topics the interested reader is addressed to the book by Michel 
( [Michel 199 1|) . In a sample of more than 300 pulsars the range of values of the surface 
magnetic field strength runs into the interval: 10.36 < logB (Gauss) < 13.33 ( [Manchester 
land Taylor 1981|) . 


Very recently, several authors ([Duncan and Thomson 1992| , [l'homson and Duncan| 
1993, [Bisnovatyi- Kogan and Mosheenko 1992|, [Bisnovatyi- Kogan 1993) have proposed two 


different physical mechanisms leading to an amplification of some initial magnetic field in a 
collapsing star. Fields as strong as B rv./ 10 14 - 10 16 G, or even more, might be generated in 
new-born neutron stars. 


According to ( [Bisnovatyi-Kogan and Mosheenko 1992| , |Bisnovatyi- Kogan 1993Q , a 
mirror-asymmetric magnetic field distribution might arise in a rapidly and differentially 
rotating proto-neutron star having, originally, both a toroidal and a poloidal component. 
The field amplification due to differential rotation leads to the formation of an additional 
toroidal field from the poloidal one by twisting of the field lines. After the first 20 seconds 
of the life of a new-born neutron star (basically, its Kelvin-Helmholtz epoch) the induced 
toroidal magnetic field could be as huge as B 10 15 - 10 17 G. 


In a second scenario (puncan and Thomson 1992| , [Thomson and Duncan 1993|) , a 


dynamo action in a differentially rotating and convective young neutron star is responsible 
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for the strengthening of some initial dipole field up to values of B ~ 10 12 — 3 x 10 13 G if 
the convective episodes arose during the main-sequence stage or to B 10 14 - 10 15 G if the 
dipole field is generated after collapse. 

In presence of a strong magnetic held the so-called neutrino-pair synchrotron radiation 
process becomes allowed : 

[e“] 4 [ e -] + i/ + P (1) 


This reaction has been studied by a few groups ([Landstreet 1967| , |Uanuto et ah 


1970|, |Yakovlev and Tschaepe 1981|, |Vidaurre 1990|, [Kaminker et ah 1991|, |Kaminker et 


al. 1992| , |Kaminkcr and Yakovlev 1993| ) for different regimes of the electron plasma. The 
calculation of the corresponding emissivity, both analytically or numerically, is far from 
obvious. In fact, this calculation appears as a multiple integral and summation over the 
variables and quantum numbers involving the wave functions of the initial and final electron 
and the corresponding statistical weights. The wave function integrals lead to Laguerre or 
Bessel functions with a complicated behaviour, so that approximations in order to simplify 
the expressions are sometimes delicate. This in fact has lead to errors in the past literature. 


In a recent paper, Kaminker et al. ([Kaminker et al. 199 1|) have studied the neutrino 
emissivity of the above process ([!]) for moderately high magnetic fields B ~ 10 12 — 10 14 G 
and high densities, in the degenerate-relativistic regime for the electrons. They claim that, 
within these conditions, the emissivity is independent of the electron density. So far, to our 
knowledge, there are no numerical tests which confirm these results. 


Given all these circumstances, together with the interest of the problem, both from 
the theoretical and astrophysical point of view, we have considered useful to reexamine 
the existing results for the above regime. Therefore, we have performed a numerical study 
of this process at high densities pY e > 10 7 g/cm 3 , where Y e is the electron fraction per 
baryon and p the matter density in e.g.s. units, and moderate temperatures T < 10 9 /l, for 
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values of the magnetic field strength B < 10 16 G. Under these conditions, the electrons are 
relativistic and degenerate. We have found a result which is in agreement with the one of 
Kaminker et al. when the electron gas is in the classical regime. We also obtain agreement 
for the corresponding analytical expressions, which we present in a simpler way than these 
authors. For the quantum regime, however, we give analytical formulae which show the 
correct dependence on B for large magnetic fields, in contrast to those given by Kaminker 
et al. This fact would be particularly important if some of the magnetic field amplification 
mechanisms described above were physically realizable in nature. 

The influence of the neutrino-pair synchrotron radiation in presence of strong magnetic 
fields merits to be examined in different astrophysical scenarios, such as the delayed 
mechanism of type II Supernovae, neutrino cooling of proto-neutron stars during the Kelvin 
epoch or the secular cooling of the neutron star. Furthermore, the neutrinos originated 
from the combined effect of one of the proposed mechanisms to enhance the magnetic field 
and the neutrino-pair synchrotron radiation process could be envisaged as a signature of 
the mechanism itself. 

This paper is organized as follows. In section 2 we discuss the calculation of the 
emissivity for the synchrotron process and present our main results. These results are 
illustrated with more detail in section 3 and 4 for two different regimes of the electron gas. 
We end in section 5 with some conclusions and remarks. 


2. CALCULATION OF THE EMISSIVITY 


The emissivity for the process (|1|) can be written as 
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£» = 3P# ES=i EKi /_ + “ dp, /,-*<„» d 3 q L, 

x Xl/(B) [1 - /(£')] 


( 2 ) 


where G is the Fermi coupling constant. The initial (final) electron, of mass m, has 
an energy E = sjm 2 + p]_ + pi (E 1 = \Jm 2 + p'l + p'£) characterized by the Landau 
quantum number n ( n') and momentum p z (p' z ) along the B-direction. We have introduced 
p\ = 2 eBn and p' 2 = 2eBn', which correspond to the classical transverse momenta of the 
electron. In the above equation, q is the four-momentum transfer and a ) = E — E' the 
energy which is carried away by the neutrino pair. f(E) is the Fermi-Dirac distribution 
function : f(E) = exp(-^^) + 1 , where fi is the electron chemical potential. The 

integration region over q is restricted by : 


2 2 
q — uj 


f > 0 . 


(3) 


The expression for A in the equation of the emissivity can be found in (|V idaurrc 


1990|, |Kaminkcr ct al. 1992|). In the relativistic limit, the relevant expression for A is : 


A = 


(Cl + c. 


EE 1 


{(-2(EE'- p z p(f 


+ (EE' - p,p()(pl - 2ql) - ^( 9 i - p\) 


T(u) 


+ Pt K EE * - pzp'z) - o (Pt ~ q!) 


$(u) 


( 4 ) 



where 
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V e " U fi( L ^-i) 2 + (^') 2 ] 

n! n 


(5) 


are Laguerre functions with argument u = ^ (qq is the component of q orthogonal 
to the magnetic field) and s = n — n'. Cy ( Ca ) is the effective vector (axial) coupling of the 
neutrino pair to the electron current, coming from both the Fierz reordered charged current 
(for electron neutrinos) and neutral current (for all neutrino species) weak interactions. 

In Eq. (f|) we have dropped a term proportional to (Cy — C\ ), which disappears in the 
extreme relativistic limit, and the interference term proportional to CvCa , which does not 
contribute to the integrated emissivity in Eq. (|2|). The argument goes as follows. The 
corresponding integrand, for CyC a, is odd under the simultaneous change of sign of the 
longitudinal momenta of the initial and final electrons, so that a symmetric integration of 
both p z and q z in Eq.(j2j) cancels this asymmetric term. We have defined p 2 = p\ + p±- If 
one takes for the electroweak mixing angle sin 2 6\y = 0.23, then Cy + C\ = 1.6748. 

One has the relationship 


u 2 -q 2 z = 2 m 2 + pi - 2 (EE' - p z p' z ) « p\ - 2 (EE' - p z p z ) ( 6 ) 

where the latter approximation corresponds to considering relativistic electrons. 

By substituting in Eq. ([|) we obtain 

A = - q\ - 9i)[p?(* - 4) - - ql - ^)«] (7) 


As mentioned above, for the range of temperatures and densities we are interested in, 
the electrons are degenerate. The product of distribution functions appearing in (|[) will 
then restrict the energies to E, E' ~ /j. Moreover, one can write the following identity : 



f(E) [1 - /(£')] = BH [HE') - f{E)] 


(8) 


where B(uS) = [exp(u;/T) — 1] _1 is a Bose-Einstein distribution function with zero chemical 
potential. From the last equation, it is apparent that the energy difference u will be 
restricted to a few times T and will be much lower than the relevant values of E and 
E', if the electrons are degenerate. Similarly, one can easily estimate that q z and q± will 
contribute as ~ T, whereas p t , p z and p' z contribute as ~ //. This allows us to neglect the 
last term in Eq. (0). In this case one gets 

A = {C \ee' A> ^ - i - "t ) 0 < 9 > 

with 0 = T — We have tested numerically that the complete expression Eq. ([4]) or 
Eq.([7|) gives approximately the same result as Eq.(|9|) for the physical conditions we are 
considering here. 

We can also perform the following approximations, in agreement with the above 
discussion : 


, _ pt (E + E’) 2 + u; 2 - 2 {pi + pf 2 ) _ (E + E’) 2 - Apl 
R U 2 eB AeB ~ 4 eB 

2 seB + p 2 z - p'? _ 2 seB + 2 p z q z 

U E + E> ~ E + E’ 1 ’ 

Due to the Pauli principle, the quantity /(E)[1 — f(E')\ will be nonzero only when E 
and E' are within a narrow interval around /i of width u ~ T. With this in mind, we have 
replaced all the slowly varying functions in (|2|) by their value around E, E’ ~ /i. 

The latter equation, together with the phase space restriction (^[), requires that q z and 


q i must lie inside the elliptical domain (|Kaminkcr ct ah 1991|) 
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( 11 ) 


where Xz = ^ and X ± = ff • 

We have numerically calculated the emissivity of the synchrotron process (|T|) for values 
of the electron density pY e > 10 7 g/cm 3 , moderate temperatures T < 10 9 K, and magnetic 
helds B < 10 16 G. The results, for this range of values, can be expressed in a compact way 
as : 


e v = 1.47 10 14 Bl 3 T 9 5 f(x) erg/cm 3 /s (12) 

where B 13 is the magnetic held in units of 10 13 G, Tg is the temperature in units of 10 9 A', 
and x = ^ is a dimensionless variable. This result is useful in order to perform analytical 
calculations. We have plotted in Fig. 1 the function f(x). As can be seen from this figure, 
f(x) first increases as x grows, and is almost constant for large values of x. This behaviour 
corresponds to different physical regimes of the electron plasma, and will be explained in 
the next sections. 


3. Large x Regime 


Let us discuss the behaviour of the emissivity corresponding to large values of the 
parameter x defined above. This situation corresponds, for example, to sufficiently high 
electron densities for a fixed magnetic held and temperature. The number of Landau levels 
which are involved in Eq. (||) is n max = On the other hand, the maximum of s = n — n' 

» 1 


can be estimated from Eq. (; 1C);) as s, 


j max ~ = x- We then have 
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for ^ >> 1. According to this idea, we have used the following approximations for the 
Laguerre polynomials in Eq. (Q) : 

* - + 

<h -> 2 J] (13) 

with J s (a ) a Bessel function and a = ^2[n + n')u. One can prove that a < s always. The 
approximation shown by Eq. © can be used if n » 1 and s « n. Because it is time 
saving and more suitable numerically than Eq. (|5|), we made use of it in our numerical 
computation of the emissivity, whenever large values of n (and s « n) were encountered. 
By changing the sum over n and n! to a sum over n' and s, and integrating the angle of q 
around the B direction, one arrives to the following expression : 


(Cf- + C 2 A ) E“,o ESi P” dp z w 

X (1 - ^ - ql - ql) B(u,) \f(E') -/(E)] 6(a) (14) 

The argument of O can be written as a = p±q±/{eB), with pj_ = \Jp 2 — pi, and c o = seB+ P~ q ~ , 
Further approximations can be made in the above equation for degenerate electrons, if one 
considers the distribution functions f(E) and f(E') as step functions in the energy. Within 
this assumption, the sum over n' can be done explicitly. One gets 

±yw-m]=£s(B + B)*!£ (i 5) 

In deriving Eq. (p~5|) , we have used the fact that s is lower than n max . By inserting the 
latter expression into Eq. ([h|) we obtain 
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e* = (Cl + C\) E“ i /_ + ~ dp z f q2 ^ 2 dq z dq ± q ± a ; 2 

x (1 - $)(u* - g - &) B(u,)Q(a) (16) 

Next we make use of recurrence relationships for the Bessel functions and obtain the 
formula : 


0(a) = 2(j;) 2 + 2( S 2 /a 2 -l)J s 2 (17) 

For the large values of s involved here, one can use the following approximation to the 
Bessel functions QGradshtcyn and Ryzhik 198C|) : 

M «)« (is) 

where z = ! and K x / 3 is the modified Bessel function, which can be further 

approximated as : 


K 1 / 3 (z) « ^exp (-z) (19) 

In this way one obtains, after some algebra, 

0(a) « -s~ 4/3 (2z) 1/3 exp (-2 z) (20) 

7r 

As can be seen from Fig. 2, the latter equation provides a reasonable approximation 
to Eq. fll7D . In this figure, we have plotted 0(a) as obtained from Eq.(^) (dotted line) and 
from Eq.(jT^) (solid line) for s = 200. Another important feature is that only values of the 
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argument a close to s will contribute. This can be understood from Eq.(^), due to the 
exponential behaviour, which effectively limits z. In fact, if we define cr as the ’width’ of 
the exponential, one has (1 — a/s) < (y) 2 ^ 3 . By substituting into Eq.([Tl]) one obtains that 
important values of q z and q± are restricted to 


( 1 _^) < (^) 2/3 

Xx S 

\q,-X*\ < A( 9l ) < V 2 XU (-) 1/3 ( 21 ) 

Pz S 

We have defined A(qq) = ^^^1 — Thus q z is restricted to a narrow interval of 
width A(qq) around Xz- This allows us to perform the integral over q z approximately. To 
the first order in A(gj_) one can write 


Xz+A(<u) 

Xz~A(q ± ) 


dq z 


2A (q A 


( 22 ) 


with the replacement q z —>■ Xz in the integrand of Eq. (|l^). This means that a ; 2 will be 
replaced by X 2 + X± ■ The integral over qq is then immediate. Since the total number of level 
differences s is large, one can substitute the sum over s by an integral over the continuous 
variable t = s/x. Thus by changing —» x / 0 °° dt and performing the remaining integrals 
one finally obtains 


e, = K (Cv + C 2 A )((5)(eB) 2 T 5 = 1.16 10 15 B 13 T 9 5 erg/cm 3 /s (23) 

7r D 

This equation implies that the emissivity does not depend on the electron density in 
this regime, in agreement with the result previously found in Ref. (|Kaminkcr ct al. 199 1[) . 


In fact, our Eq. (|23|) is close to the one obtained in this reference. It is also in good 
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agreement with the values of f(x) obtained numerically (and plotted in Fig. 1), as can be 
seen by comparing Eq.(p3|) with Eq.(|T^). 


A numerical fit which reproduces the behaviour of f(x) for x > 2 to better than 4% is 
given by 


-5.0224 - 8.1289a; + 9.2892a: 2 
1.0293 + 2.0605a; + 1.0727a; 2 


(24) 


4. Low x Regime 

We now address the question of whether an increase of the magnetic held will always 
give a larger neutrino emission. In Fig. 3 we present the corresponding emissivity (in cgs 
units and logarithmic scale) for a fixed temperature Tg = 1 and four values of the electron 
density pY e = 10 9 ,10 11 ,10 13 , 10 u g/cm 3 , as a function of Bi 3 (R 13 ranging from unity up 
to 10 3 ), as obtained numerically. For these high magnetic fields, the number of populated 
Landau levels (n max , defined above) can be of order unity, and we enter into the quantum 
regime. We have used, in these cases, the expression of A as given by Eqs. (fD and (||) 
directly, instead of making the approximations shown in section 3. 

As can be seen from Fig. 3, for a given density and temperature, the emissivity first 
increases and, after reaching a maximum value, will fall to zero for large values of the 
magnetic held. This can be understood since the number of possible n —> n' transitions 
decreases as ^ S> 1. In fact, a rough analytic expression in this region can be obtained by 
putting s = 1 in Eq. ( ]T4|) and u — Therefore one has : 


T ~ 1 


rv 


0 


(25) 
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The integrals in Eq. (0 can be done analytically and one obtains the following 
expression : 


e v = 6.6 10 12 £q 3 Tg x 5 e l ^ x erg/cm 3 / s 


(26) 


We have verified that the latter expression gives values which are in agreement with 
our numerical results around the maximum of the cmissivity. This is shown in Fig. 4, 
where we have compared our results for pY e = 10 11 g/cm 3 (solid line) with the prediction 
of Eq. ( |26| ) (dotted line). We have also plotted (dashed line) the corresponding analytical 
approximation given by ([Kaminker et al. 1991|) for this case. As can be seen from this 
figure, Eq. (|26]) provides a reasonable approximation to the emissivity, which works better 
than the formula of Kaminker et al. 


5. Comparison with other processes 


In order to investigate the relevance of the process studied here for neutron star cooling, 
we have made the comparison with the emissivities corresponding to other neutrino processes 
which are competitive with the synchrotron emission. We have considered pair production 
e + e _ —> vu ,plasmon decay T —> vv, bremmstrahlung e~(Z,A) —> e~{Z ) A)uu and 
photoproduction ye - —» e~vv. For these processes, we have made the assumption that 
they do not vary significantly with the magnetic field. The numerical fit to these emissivities 
have been taken from [Munakata et al. 1985 . Although these rather simple formulae are 
not the most up-to-date available fits to the above processes, they serve to our purpose as 
a first approximation (see, for example, [Itoh et al. 1989|, for a more elaborated fit). The 
bremmstrahlung process is taken from [Maxwell 1979| . A more complete calculation, as in 
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Itoh et al. ( itoh ct al. 1989 ), gives the same order of magnitude in the region where this 
process dominates. 

In Fig. 5 we present the result of comparing all these processes for a temperature 
T = 10 8 /i as the product pY e varies from zero up to 10 12 g/cm 3 . The bremmstrahlung 
energy emission was calculated assuming that the dominant nucleus is 56 Fe. The 
synchrotron emissivity is plotted (solid lines) for two values of the magnetic held : B 13 = 10 
and h>i 3 = 100. Other relevant processes are : bremmstrahlung (dashed-dotted line), 
plasma (long dash), and photoneutrino emission (short dash). As can be seen from this 
figure, the synchrotron emission is competitive with the above processes within the range 
pY e 10 9 — 10 12 . Moreover, as pointed out by Pethick and Thorsson QPcthick and Thorsson| 


19931) , band-structure effects can suppress bremmstrahlung by a factor of 10 or more for 


temperatures less than about 10 9 /l. In this case, the synchrotron emission would be the 
dominant process in the above electron density range, if the magnetic held reaches values 
of the order ~ 10 14 G. 

For higher temperatures, the synchrotron emission corresponding to a given value of 


B ’switches on’ at lower electron densities, as can be inferred from Eq. (p 6 |). However, 
other processes have a faster increase with temperature and, therefore, the dominance 
of synchrotron emission reduces to a narrow interval of densities, although it effectively 
competes for high densities. This is shown in Fig. 6 , where we have made the above 
comparison for a temperature T = 10 9 K (pair emission is represented by the dotted line). 


6. CONCLUSIONS 

We have performed numerical calculations of the synchrotron emissivity from 
relativistic degenerate electrons. This calculations allow us to present the results in terms 
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of an universal function f(x) which can be used in astrophysical codes. For two different 
regimes of the electron gas we have derived analytical formulae, in a simpler way than 
previous references. These formulae have been tested, and we have found a reasonable 
agreement with our numerical calculations. We also have compared our results to the 
analytical formulae derived recently by Kaminker et al. Agreement is found for the classical 
region, where many Landau levels contribute to the emissivity , but some differences arise 
in the quantum regime, for the values of the magnetic field recently suggested in new-born 
neutron stars. 

We have shown that neutrino-pair synchrotron radiation for moderately high magnetic 
fields B > 10 14 G is an efficient cooling mechanism for temperatures not larger than 10 9 /l, 
and can compete effectively (or even dominate) with other processes. 

We claim that the influence of the neutrino-pair synchrotron radiation in presence of 
strong magnetic fields (~ 10 15 G) merits to be examined in different astrophysical scenarios, 
such as the delayed mechanism of type II Supernovae, neutrino cooling of proto-neutron 
stars during the Kelvin epoch or the secular cooling of the neutron star. Furthermore, these 
neutrinos originated from the combined effect of the dynamo action with the neutrino-pair 
synchrotron radiation process could be envisaged as a signature of this mechanism. This 
will be the subject of future investigations. 
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Figure Captions 

Figure 1.- The function f{x) appearing in Eq. flT2|). See the text for the definition of 
the dimensionless variable x and comments about its behaviour. 


Figure 2.- The function 0(a) as obtained from the approximation Eq.(|20j) (dotted 
line) compared to Eq.(]T7|) (solid line) for s = 200. 

Figure 3.- Neutrino synchrotron emissivity as a function of the magnetic field B 13 for 
different values of the electron density. The temperature is the same (T = 10 9 /i) in all 
cases. 

Figure 4.- Comparison of our numerical results for pY e = 10 11 g/ civ? (solid line) with 
the analytical approximation Eq. (^6|) (dotted line). We have also plotted (dashed line) the 
approximation given by ( |Kaminker et ah 1991|) for this case. 

Figure 5.- Competition of synchrotron neutrino emission (solid lines) with other 
processes, as a function of the electron density, for a temperature T = 10 8 /l. Two 
different values of the magnetic field (.B 13 = 10 and Bi 3 = 100 have been considered. The 
bremmstrahlung emissivity (dash-dotted line) has been calculated for 56 Fe. Plasma process 
is represented by long-dashed line, and photoneutrino by short-dashed line. 

Figure 6.- Same as figure 5 for T = 10 9 /i. Dotted line corresponds to pair emission. 
















